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Abstract 

We show that a class of spin models, containing the Ashkin-Teller model, 
admits a generalized random-cluster (GRC) representation. Moreover we 
show that basic properties of the usual representation, such as FKG inequali- 
ties and comparison inequalities, still hold for this generalized random-cluster 
model. Some elementary consequences are given. We also consider the duality 
transformations in the spin representation and in the GRC model and show 
that they commute. 
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The introduction by Fortuin and Kasteleyn ||FK| , pl| |F2|| of the random-cluster 
model in the late 60s has given rise to numerous important results. First it provided 
a unified representation of several famous models, including the Ising, Potts and 
percolation models, thus allowing the comparison between them. It also brought a 
whole class of models interpolating between the latter ones. The random-cluster 
representation has been used in many recent proofs in statistical mechanics, for ex- 
ample in large deviations theory Q . The fact is that this model has several nice 
properties, as FKG and comparison inequalities, allowing to derive non-perturbative 
results for the original models. One of the properties which has also often been used 
is that the two-dimensional random-cluster model is self-dual, and that this duality 
commutes with the duality of the original models; this has been used for example 



in the study of the decay of the connectivity in the Ising model ||CCS|| . Other ap- 



plications of this representation have been found in numerical studies, in particular 
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the Swendsen-Wang algorithm is based on it. 

It would then be interesting to be able to extend this representation to a wider class 
of models, while keeping most of its properties. This appears to be possible. We 
show that the Ashkin-Teller model (and a class of models generalizing the Ashkin- 
Teller model, and containing the partially-symmetric Potts models) admits a similar 
representation, which in fact generalizes the usual one. The nice point is that it is 
still possible to prove FKG inequahties, comparison inequalities and commutativity 
of the dual transformations for this new representation. 

Such a representation has already been considered in |[WD| , |SS|| . The main goal in 



these papers is to develop a Swendsen-Wang type algorithm for the Ashkin-Teller 
model. A closely related representation has also appeared in the study of partially 



symmetric Potts models ||LMaR|| . Their representation appears as a special case of 
the one studied here. Nevertheless, properties of the measure were not studied in 
these papers. 

Although the Ashkin-Teller model has been introduced more than half a century 
ago [|AT|| , there are still several open questions about this model. Some of the tools 



developed for the study of the Potts model via the random-cluster representation 
are useful in the study of the Ashkin-Teller model. In this paper, we focus on the 
properties of the two-dimensional model, and give only some elementary applica- 
tions of the inequalities. At the end of the paper, we discuss possible extensions 
of the results. We shall consider more elaborate applications in a separate publi- 
cation. One of the main points of the paper is to show that elementary methods 
can be used to study the duality transformation of the spin model and the random- 
cluster representation. It is advantageous to derive the duality transformation using 
the high-temperature expansion based on the elementary formula ( |2.5|) ; moreover, 
this approach allows to study correlation functions and boundary conditions very 
explicitly. The random-cluster representation is not more difficult than the high- 
temperature expansion; it is based on the elementary formula ( [3.221 ). 



After we finished this work we received the paper ||UM|| by L. Chayes and J. Machta. 
In this paper graphical representations are developed for a variety of spin-systems 
including the Ashkin-Teller model. These representations are used in connection 
with Swendsen-Wang type algorithms. The case of the Ashkin-Teller model is stud- 
ied in details. Although the presentation of the model is different (compare e.g. 
the phase diagrams), essentially all our results about the random-cluster model are 
explicitly derived in | ](J1V1| | (see in particular Propositions 3.5 and 3.6 therein). 

Acknowledgements: We thank L. Chayes for discussions and communicating us 
his results with Machta. We also acknowledge discussions with L.Laanait and J.Ruiz 
about the duality transformation. 



1 The Ashkin-Teller model 



Lattices and cell-complexes 

The model is defined on 2^ or on some bounded subset A C Z^, 



(1.1) 
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We call sites t the elements of the lattice 2^. Two sites t and t' are nearest- 
neighbours if Iti— t']^| + |i(:2— ^2! — 1- By definition the boundary of a site is the empty 
set. We call bonds b =<t,t'> the subsets of IR^ which are straight line segments 
with the nearest-neighbours sites t and t' as endpoints. The boundary of a bond is 
5b = {t, t'}, and the boundary of a set B of bonds is the set 5B = {t E : t E 5b 
for an odd number of bonds b E B}. Finally we call plaquettes p the subsets of IR^ 
which are unit squares whose corners are sites. Their boundary is the set of the four 
bonds forming their boundary as a subset of H^. With this structure the lattice 
becomes a cell-complex, which we denote by IL. 

Another lattice is important, the dual lattice (Z^)*, 

(Z^)* := {t = {tiM) : + 1/2 G Z, z = 1, 2} . (1.2) 

We can of course define the same objects as before for the dual lattice, they will be 
denoted t* , b* and p* respectively. The dual cell-complex will be denoted by IL*. 
The following important geometrical relations hold: 

1. each site t is the center of a unique plaquette p*, 

2. each bond b is crossed by a unique bond b*, 

3. each plaquette p has a unique site t* at its center. 

A subset A C Z^ is simply connected if the subset of IR^ which is the union of 
all plaquettes p*{t), t G A, is a simply connected set in IR^. 

Dual of a set 

Let A C Z^; we will also denote by A the following subset of IL: the sites of A are the 
elements of A (as subset of Z^); the bonds of A are the bonds of IL whose boundary 
belongs to A; the plaquettes of A are the plaquettes p of IL whose boundary is given 
by four bonds of A. We will denote by B{A) the set of bonds of A. 
We now define a dual set for A. We will define another notion of dual set later (see 
subsection p.3|) . 

We define A* C IL* in the following way: the plaquettes of A* are all plaquettes of 
IL* whose center is some site of A; the bonds of A* are all bonds of IL* belonging to 
the boundary of some plaquette in A*; the sites of A* are all sites of IL* belonging 
to the boundary of some bonds in A*. 

Configurations, Hamiltonian and Gibbs states 

A configuration uj of the model is an element of the product space 

n:=[{-l,l}x{-l,l}f\ (1.3) 

The value of the configuration to = (cr, r) at t G is uj(t) = (cr(t), T(t)). 

Let A C Z . A configuration cu is said to satisfy the (+, +)-boundary condition 

in A if 

w(t) = (l,l) Vt^A. (1.4) 
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The Ashkin-Teller Hamiltonian on A is 



Ha = - {Jo(^i(^j + JrTiTj + J^r<yiajTiTj} , (1.5) 



<i.j>: 

{ij}nA7^0 



where J^; Jr and J^j. are real numbers called coupling constants. 

The Gibbs measure on A with (+,+)-boundary condition is the probability 
measure given by the formula 



S(+.+)(A)"iexp(-iJA(cj)) if uj{t) satisfies the (+, +)-b.c. on A, 

otherwise. 

(1-6) 

where the normalization S^"'"'+)(A) is called the partition function with (+,+)- 
boundary condition. 

In the same way, we can introduce (+, — )-, (— , +)- and (— , — )-boundary condi- 
tions by imposing the corresponding value to u outside A. 
Notice that the Ashkin-Teller model has the following symmetries : 

/^A^((t^> t)) = /^a"((^> -^)) = /^A^((-^> -r)) = /^a"((-^, -r)) , (1.7) 
so we consider only (+, +)-boundary condition. 

We also define the Gibbs measure on A with free boundary condition 

/i((u;) := S^(A)"^ ]J exp{ J^aiCTj + Jr^Tj + JarCTiajTiTj} , (1.8) 

<i,jr>CA 

where the normalization ^•'^(A) is called the partition function with free bound- 
ary condition. 



Remark: For J^^- = 0, the Ashkin-Teller model reduces to 2 independent Ising 
models, while for J„ = = it becomes the 4-states Potts model. 



We will always suppose that the coupling constants satisfy 

Ja>Jr>Jar- (1-9) 

Note that there is no loss of generality in doing this choice. Indeed we can always 
transform ( |1.5|) to obtain this order. For example, if Jg-r > Jr, then we can make 
the following change of variables: (cxj, Ti) (cTj, 9i), where 9i = aiTi. 

In this paper, we further impose that 

Ja > 0, Jr > 0, tanh Jo-T > — tanhJ^-tanhJ,- . (1-10) 
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2 Duality of the Ashkin— Teller model 

Duality of the Ashkin-Teller model has been known for a long time ^ . However, 
for the sake of completeness, as well as to fix the notations which will be used when 
considering the duality of the random-cluster model, we give here a straightforward 
account of this transformation. 

2.1 Low temperature expansion for (+,+) -boundary condi- 
tions 

Let A C be bounded and simply connected. Let us now consider the Ashkin-Teller 
model defined on A, with (+, +)-boundary conditions and with coupling constants 

With this kind of boundary conditions, we can describe geometrically all configura- 
tions (cr, r) of the model by giving the sets 

M^:={p*{t):teA,at = -l}, 
Mr:={p*it):teA,Tt = -l}. 

The boundaries of these sets, considered as subsets of IR^, define two sets of bonds 
of IL*. Maximal connected components 7^, 7r of these sets of bonds are called a- 
and r-contours respectively. We will call closed contours contours such that 
^7 = 0. The length of a contour is its cardinality as a set of bonds and is denoted 
by I7I. A configuration of contours is a set of closed contours such that: (a) any 
two cr-contours are disjoint (as sets of bonds and sites); (b) any two r-contours are 
disjoint (in the same sense). (There is no constraint between the a- and r-contours.) 
Such a set will be denoted by (7^,7^), where 7^ denotes the set of a-contours and 
7^ the set of r-contours. To each spin configuration u = (o", r), it is possible to 
associate a unique configuration (7^, 7^) of contours. 

Remark: If A is simply connected, then the converse is also true. If it is not simply 
connected, then it will generally be false. Indeed, suppose A is a square with some 
hole in it, with (+, +)-boundary condition. Then only configurations of contours 
such that there is an even number of a (and r) -contours winding around the hole 
correspond to some spin configurations. This will be important when considering 
duality. 

Let us now introduce the weights of contours 

l^aila) ■■= exp(-2(j^ + JaT)\la\) , ^^(7^) := U^^^X^^Al^) , 
^r(7r) := exp(-2(j^ + Jar)\lr\) , ^Al^) ■= U^^ex^ ^rilr) ■ (2.1) 

Introducing the following interaction between the contours, 

UJarila, It) ■= eXp(4J^^|7^ f] 7^ | ) , (2.2) 

where [7^ H 7^| is the cardinality of the set of bonds belonging simultaneously to 7^ 
and 7^, the partition function in A with (+, +)-boundary condition can be written 

S^+'+^A = E ^.(7j^.(7j^..(7,,7j , (2.3) 

— <T — T 
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where Ci is some constant depending on A but not on the configurations which does 
not affect the results below. The sum is over families of closed a- and r-contours. 

Remark: If J^^- > the interaction is such that the a- and r-contours will attract 
each other while they will repel each other when J^-r < 0. 

It is therefore natural to use a normalized partition function with (+, +)-boundary 
condition which is defined as 

:= ^ OO^ij^P^xMril^^lr) ' (2-4) 



2.2 High-temperature expansion for free boundary condi- 
tions 

Suppose A C 2^ is bounded and simply connected. Let A* be the dual of A as de- 
fined earlier. We consider the Ashkin-Teller Hamiltonian on A* with free boundary 
condition and coupling constants J*, J* and J*^. 

We now proceed in doing a high-temperature expansion of S-^ 

S'^ = ^ ]^ (coshJ* + crjCrjSinhJ*)(coshJ* + rjTjSinhJ*) X 

o-,r <iJ>cA* 

X (cosh J*^ + aiCTjTiTjSmhJ*^) 
= (coshj;coshj;coshj;j 1^(^)1 1111(1 + aicx^tanhj;) x 

o-.T <i,j> 

X (1 + TiTjt8iBhJ*){l + aiajTiTjianhJ*^) . (2.5) 

Defining 

s = tanhj;, t = tanhj;, / = tanhJ^^ , (2.6) 

and 

_ s + tl _ t + sl _ I + St 

^ - 1 + str 1 + str 1 + str ^ ' 

the above sum becomes 

(1 + st/)l^(^'l E n {1 + -^^i^J- + TnTj + La.ajTiTj} . (2.8) 

0",T" <i,j> 

Expanding the product, we obtain a sum of terms that can be indexed by {ri^,T]T-), 
Va-, Vt £ {0, l}^*^^) (we recall that B{A) is the set of bonds of the cell-complex A). 
This is done in the following way: 



1. Each time we take one term 1 in ( |2.8| ), we set 

??a(<^,J>) = 0, r]r{<ij>) = 0, 



2. Each time we take one term Sa^aj in (|2.8|) , we set 



Va{<i,j>) = 0, r]r{<i,j>) = 1, 
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3. Each time we take one term Taiaj in ( p^.8D , we set 

ri^{<i,j>) = 1, r]r{<i,j>) = 0, 



4. Each time we take one term LaiCjjTiTj in (|2.8|) , we set 



Va{<i,J>) = 1, Vr{<i:j>) = 1- 

To each of these pairs (r/o-,-?]^) we associate a configuration of a- and r-contours 
(7^,7^), where the ■jo- are maximal connected components of {b G B{A) : rj^j^b) = 1} 
and 7t- are maximal connected components of {b G i3(A) : rir{b) = 1}. 



Note that we have interchanged a and r, for later convenience, see section pT3 . 
We now sum over a, r. Using the fact that o",^*'^^ = Er t"?''"^^ = 0, VA; G IN, 
we see that the only contributing configurations are those with only closed a- and 
r-contours. We obtain 



^l + St/)|e(A)|4|A| J2 {S^lr\~\l.^lr\T\l.\-\la^lr\L\l.^lr\y (2.9) 



where I7I denotes the cardinal of 7, considered as a set of bonds. We define the 
normalized partition function with free boundary conditions to be 

Z^A*:= i^S^lJ-K'^Vr^l.^-K'^VL^Zr'^lr^^ . (2.10) 



2.3 Duality 

Proposition 2.1 Let A be a simply connected bounded subset of Z^. Let V = 
{{x,y,z) G '■ X > y > z,y > 0,tanh2; > — tanhxtanhy}. Let {Jo-, Jr, Jar) E V be 
the coupling constants of the Ashkin-Teller model defined on A with (+, +)-boundary 
conditions. Then the following relations 

5(J:,J;,J:J = exp(-2(J, + J..)), 

t{j:,j:,j:^) = exp(-2(j. + j,.)), (2.11) 
= exp(-2(j. + j.)) , 



(where S, T and L have been introduced in ( \2. define a bisection from T> on itself, 
such that 

Z^^AiJa, Jr, Jar) = A^J*., J.*, Jtr) ■ (2-12) 

On the closure ofV, the application is still well-defined, but takes values in IR^ and 
is no more everywhere invertible. 



The proof is straightforward algebra; it is given in the appendix. 
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2.4 The self— dual manifold 

Proposition 2.2 The self-dual manifold, i.e. the set of fixed points of the duality 



relations ( \2.1]\ ), is given by 



l= ]-''-'-\ (2.13) 
1-st + s + t ^ ^ 

where s = tanhJ^-; t = tanhJ,- and I = tanhJo-r- 

Proof. We want to find the values of J„, Jr and Jo-r such that J* = J^, J* = Jt 
and J*^ = J„r- In particular one must have 

1 + stl + + ^ ' 

We have used 

g-2(x+j/) ^ (1 -tanhx)(l -tanhy) 
(1 + tanha;)(l + tanhy) 

After some algebraic manipulations, (|2.14|) can be seen to be equivalent to 

1= ^-''-'-' . (2.16) 
1-st + s + t ^ ^ 

The two other relations are seen to be satisfied for these values of / by substitution, 

s + tl _ {l-t){s + t) _(l-t)(l-/) 
1 + stl ~ (l + t)(l-st) ~ (1 + t)(l + /) ' 

t + sl _ (l-s)(s + t) _(l-s)(l-/) 
l + stl ~ (l + s)(l-st) ~ (l + s)(l + /) ■ 

□ 

Remark: Note that, in contrast to the 2 dimensional Ising model, this self-dual 
manifold does not coincide with the critical manifold [ |W| , pi[ . For example, in the 



= plane, the self-dual line and the critical line coincide only when Jq.,- < J^, 
then the critical fine sphts into 2 dual components. See section (|4.2| ) for an estimate 
on the location of these lines. 



3 The Random— cluster model 

In this section, we introduce the generalized random-cluster model (GRC) and show 
its connection to the usual random-cluster model and to the Ashkin-Teller model. 

We introduce the model by discussing successively the configuration space, the a 
priori measure (generalized percolation measure) and the generalized random-cluster 



measure. 
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3.1 The model 

Configuration space 

For every bond b, let := {0, 1} x {0, 1}. 

The configuration space is the product space T := Tf^^ \ where is the 

set of bonds of Z^. 

A configuration of bonds n is an element of the configuration space. The 
value of the configuration n = {11^,11^.) at a bond b will be denoted either n{b) = 
{n^{b),nr{b)), or Ub = {n^,b,nr,b)- 

Bonds b such that na{b) = 1 are said to be cr-open, while bonds b such that na{b) = 
are said to be u-closed. In the same way we define r-open and r-closed bonds. 

If n = {n^,n^) is some configuration of bonds then n is the configuration given by 

nb = {1 - n„M, 1 - nr,b)- 

Let n G T. We define a notion of connectedness for sites, given the configuration n. 
The site i is cx-connected to the site j, given the configuration n, if there exists 
a sequence to = i,ti, ...,tk-i,tk = j of sites such that n(^(< ti,ti^i >) = — 
0,...,A;-1. 

Maximal connected components of sites are called cr-clusters. The number of a- 
clusters in a configuration n which intersect a given set A is denoted by N^{n\A); 
note that each isolated site is a cluster. 

Two sets are cr-connected, given a configuration n, if there is a point of the first set 
which is (T-connected to a point of the second set. 
If i and j are cr-connected, we will write 

. <T 

We make the corresponding definitions for r. 

The a priori measure 

On the configuration space we introduce an a priori measure, which we call gener- 
alized percolation measure (GP measure). 

We introduce for each b e i3(Z^) a probability measure A;, on T;,, given by 

A;,((0,0))=ao(6), A5((l,l)) = a,,(6), 

\b{{l, 0)) = a,{b) , A5((0, 1)) = ar{b) . (3.2) 

Let i3 be a finite subset of B{W?). The generalized percolation measure in B 

is defined as the following product measure on T 

XB{n)= n «o(&) n n «r(&) n «-(&). (3.3) 

beB: beB: beB: beB: 

ni,=(0,0) "(, = (1,0) nb=(0,l) "(,=(1,1) 

The generalized random-cluster measure 



(3.1) 



Random-cluster representation of Ashkin-Teller model 



10 



Let A be a bounded simply connected subset of 2^. 
We recall that 

B{A) -.^{beJL-.SbcA}. (3.4) 
We introduce another set of edges associated to the set of sites A, 

B+(A) := {6 e IL : 56nA7^ 0}. (3.5) 

We introduce two kinds of boundary conditions. 

The configuration n satisfies the (+, +)-boundairy condition on A if 

n{b) = (1,1) V6^i3+(A). (3.6) 
The configuration n satisfies the free boundeiry condition on A if 

n{b) = (0,0) V6^i3(A). (3.7) 
Notice the fact that the set of bonds in each of these definition is different. 

Remark: 1) (+, +)-boundary condition corresponds to what is usually called 
wired boundary condition. 

2) We can also define more complicated kinds of boundary conditions 
by imposing the corresponding values for the configuration outside B(A) or B^{A). 

We introduce the following notations 

n:n satisfies the (+,+)-b.c. on A 
S/,A • satisfies the free b.c. on A (^•^) 

The generalized random— cluster measure with (+, +)-boundary condition 
on A is the probability measure on T given by 

^ "^^'^^ ~ ,\ N^(n \ l) iv^(» | A) if n satisfies the (+, +)-b.c. on A, 
otherwise. 

(3.9) 

where and Qt are two positive real numbers. 

The generalized random— cluster measure with free boundciry condition on 

A is the probability measure on T given by 

Er^U)(>'- ' %^(z.|A) if^^-tisfiesthefreeb.conA, (3,0) 
otherwise. 

Relation to the usual percolation and random-cluster measure 

For special classes of functions, which we define below, the expectation value in the 
GP or GRC measures can be related to the expectation value in some percolation 
or random-cluster measures. 
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We introduce three classes of functions on T. 

Let J-'^ be the set of functions on {0, 1}^^^ \ and JF^ be the set of functions on 
[{0,1} X {0,1}]^(^'). We define 

:= {/ G ^2 . ^ jri ^ith f{n) = U{n^) Vn} , 

J^r := {/ G : 3/. G J^' with /(n) = /.(nj Vn} , (3.11) 
:= {/ G ^2 . gj^ ^ ^1 ^itj^ y(^) ^ ^^(^^ V n,) Vn} . 

where (rz^ V?i^)(6) := max(no-(6), ^,-(6)). 

We prove now an elementary lemma relating the GP measure on B to the usual 
percolation measure on B, which is defined on {0, l}^*^^ ^ by 

CB{n\p):= llPib) n(l-PW)' (3-12) 

beB: beB: 
n^ — 1 

where < p{b) < 1 V6. 

Lemma 3.1 1. If f E then C,B{fa\v'^) = «<t(&) + a<7r(&)) = >^B{f)- 

2.IffE Tr then Cjs^frW) = a^b) + a^Ab)) = Ab(/). 

5. /// G t/ien (isifblpib) = a^b) + a^b) + a.,(6)) = Ae(/). 

Proof. We have (omitting the dependence on b of the probabilities) 



"r fees 



fee/3 n^,b=±l 



cr,b 
(TT 



beB 



E /^(i^a) n ("-^ + ^'^^) n (^0 + ctr) 



66B: 66B: 



□ 



The two generalized random-cluster measures are also related to the corresponding 
usual random-cluster measure on A, which are defined on {0, 1}^(^') by usmg 
notations similar to ( p.8|) ) 

I otherwise. 

p{{n\p, q) := ( e'TSIiS^^'^ ^ ^ ^^'^ ^' (3.13) 

1 otherwise. 

where A^(n|A) is the number of clusters in n intersecting A. This is proved in the 
following 
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Lemma 3.2 

/ G JF^ i^liflqa, 1) = plifalp = a<T + a<^r, Qa) , 

where o means free (resp. wired) boundary condition for the usual random-cluster 
model, and • means free (resp (+, +) ) boundary condition for the GRC model. 



Proof. As G J-'a we have, by lemma (and omitting the dependence on b of 
the probabihties) 



Eo,A/(ii)AB°(A)(n)g^ 



Na{n\A) 



Eo,A Ae°(A) (njqa 

E.,A faill)CB'(A)ill.\P = 0-0 + aaT)qa 

CB'{A){n\p = a^ + a^r)q^^^^^ 

PliUlP = «<7 + aar^qa) ■ 



N{n\A) 



□ 



3.2 Relation to the Ashkin— Teller model 

The Ashkin-Teller model defined in section || and the generalized random-cluster 



model defined in the section ^TT] are closely related as is shown in the following 
Proposition 3.1 Let 

a„ = e"^-^"(e"^'^"" - e"^-^") , 
Qr = e"^-^"(e"^'^"" - e"^"^") , 

a^r = 1 - e^2(^'^+^-) - e-2(^^+-'''^) + e-2(-^-+^-) . (3.14) 



The constants define a probability measure ( ^. ^ ) on {0, 1} x {0, 1} 

if, and only if 

Jq- ^ J(JT ^ Jt ^ Jar 1 

Jo- > , Jr > , tanhJo-T- > — tanhJ^-tanhJi- . 

Moreover, with this choice of probabilities, 
1. 

Z++ = C2 E+,A Ae+(A)(n)2^'^(^l'')2^^(^l'^) , (3.15) 

Zf = CsEm AB(A)(n)2^-(^l^)2^^(i^l^) • (3-16) 



/^r(fXArB) = J/^(«;>^|2,2), (3.17) 
/^((^ArB) = ^^{(«:M|2,2). (3.18) 
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where C2, C3 are some constants independent on the configuration and k'^ is the 
characteristic function which is one on the configurations such that no finite a- 
cluster contains an odd number of sites of A; n\ is the corresponding characteristic 
function for r. Finally, a a '■= DigA^i '^'^'^ '^b '■= Hies'^i- 



Proof. 1. Let us first show that define a probabihty measure on 

{0,l}x{0,l}. 

By definition, oq + Oo- + + a^r = 1- Hence we just have to check their positivity 
Evidently ao > 0. 

> > Jar , 

Qr > Jr > Jar , 

and 

_2j 1 + e"^('^'^+'^^) 1 - e-2-^^^ 1 - e-^-^^ 1 - e'^-^^ 
a^^>0^e < ^ ^-2J. ^ ^ ^ g-2j.. > - ^ ^ g-2j. ^ ^ ^-2J. ' (3-19) 

but this is just tanhJo-r ^ — tanhJo-tanhJ^. 
Now note that 

JcT > Jar , Jt > Jar , tauhj^-r > — tauh J^-tanh ^ ^(T > , Jr > . (3.20) 

Indeed, suppose Jcr < and Jt > 0, then Jg-r niust be positive and therefore larger 
than Jo- which is a contradiction. If J^- < and < then in this case 

tanhlJo-rl < tanh| Jo.|taiih|J.r| < (tanh| Jq-tI)^ , (3.21) 
which is also a contradiction. 

We show that the partition functions can be expressed in term of the denominator 
of the generalized random-cluster measures. 

The weight in the partition function can be expanded in the following way 
exp{J^aiaj + Jr^Tj + J^TO-iCJjTiTj} 

= CeXp{(Jo + J^T)i<yiCrj - 1) + {Jt + JaT){riTj - 1) + Jario-iCrj - l){TiTj - 1)} 
= C{ao + a„5„^„^ + aJr^r^ + a^r^a.a/r.rj } , (3.22) 

where C is some constant. The partition function with free boundary condition can 
then be written 



f,A b=<i,j>eB{A): 6=<i,j>G8{A): 

"it, 6 = 1 "T,b = l 



= C3^Ae(A)(R)2'^-(^l'')2^^(^l^). (3.23) 

The case of (+, +)-boundary condition is treated in the same way. 
2. We finally prove ( CT) and ( CT )- 

The same expansion as above can be done on the correlation functions. We then 
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E/,A AB(A)(ii)2^-(^il^)2^-(iil^) 
E/,AAg(A)(n)/^^(ll)/^^(ll)2^-(^l^)2^-(^l^) 

E/,AAB(A)(li)2^-(^^l^)2^-(lil^) 

z/{(/s:>^|2,2), (3.24) 

where we have used the fact that the only configurations that will give a non zero 
contribution must be such that cr^ = = 1, Vcr, r. But this is only possible if the 
intersection of A and any cr-cluster, as well as the intersection of B and any r-cluster 
contains an even number of sites. 

The case of (+, +)-boundary condition is treated in the same way, using also the 
fact that the sites belonging to the infinite cluster have the fixed value (1,1). 

□ 

Remarks: 1. As already stated previously, if the conditions on the order of the cou- 
pling constants of the Ashkin-Teller model are not satisfied, then it is still possible 
to use the random-cluster representation. Indeed, by first doing a change of vari- 
ables, we can always write down the model in the required form. As an important 
example, consider the case = Jr ^ Jar- The change of variables (a, r) i-^ (6*, r) 
results in Je > Jr = Jqt to which we can apply the GRC representation. Notice that 
the resulting random-cluster measure has the property that a,- = 0. This implies 
that the ^-open bonds play the role of the (random) graph on which the r-bonds 
"live" . We will return to this particular case later. 

As an important particular case of the above proposition, we have 

^^l^{cr^) = ly+{^^A''), (3.25) 
^t+{a,a,) = ut{t^j). (3.26) 



obtain 



2. In fact the generalized random-cluster model with q^, Qt can also be related 
to some spin model. More precisely, if we consider the model defined in the following 
way: 



n 



(3.27) 



where cTj G {1, go-} and Tj G {1, ^r}, then an analogous proposition as the one 
above still holds. These models are usually called (go-, Qt) -cubic models [?]; they 
may be thought of as resulting from two coupled Potts models. In the case Jr = Jar 
we recover the partially symmetric Potts model |DLMMR, LMaR |. Notice also that 



the Hamiltonian ( |3.27| ) cannot be cast into the form of the Potts models considered 
by Grimmett Q. 

3. More complicated boundary condition can be treated in exactly the same way as 
here. 



We are now going to show that the generalized random-cluster model is self-dual. 
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3.3 Some geometrical results 

In section |I| a definition of dual set was introduced. Tfiere, tfie relevant variables were 
spins and so the primary geometrical objects were sites. The dual set was therefore 
defined starting from those sites, building the corresponding dual plaquettes and 
completing the cell-complex thus obtained. 

We are now going to give another notion of dual of a set. As in the random-cluster 
model the variables are the bonds, it will be natural to build the dual set starting 
from the dual objects associated with the bonds, namely the dual bonds. 

Let i3 be a set of bonds. We define an associated cell-complex A(i3): its set of bonds 
Ai{B) is the set of the bonds in B; its set of sites Ao(-B) is the set of the boundaries of 
these bonds; the set of plaquettes A2{B) is the set of the plaquettes whose boundary 
belongs to B (there may be none). 

Let B be such that Aq{B) is a bounded and simply connected. 
We define the dual of the set of bonds B: 

B = {bel.* -.b crosses some b e B} . (3.28) 

The corresponding dual cell-complex is A := A{B). 

Let n G {0, 1}^ be a configuration of bonds. We define the dual configuration 
n e {0, 1}^ to be 

ng = l-nfe. (3.29) 
where b is the bond of B intersecting b. 

For a given configuration of bonds n we denote by (A, n) the graph whose vertex 
are the sites of A{B) and whose edges are the open bonds of n. 

We then have the following two relations: 

N{n) = \A\ - \n\ + l{n) , (3.30) 
N{n) = l{n) + 1 , (3.31) 

where N{n), |A|, |n| and l{n) are respectively the number of connected components, 
the number of vertices, the number of edges and the cyclomatic number]^ of the 
graph (A, n); /(n) is the cyclomatic number of the graph (A,n). 

Relation (|3.3CI| ) is just the well-known Euler formula for the graph {A,n) and can 
be easily proved (see, for example. Theorem 1 in ||Be|| ). Relation ( |3.31| ) becomes 

^An elementary cycle of an oriented graph (V, E) (i.e. a graph whose edges have an orien- 
tation) is a sequence of distinct edges (ei, . . . , e„) such that every ek is connected to ek~i by one 
of its extremities and to Cfc+i by the other one ( eo := e„, Cn+i := ei) and no vertex of the graph 
belongs to more than two of the edges of the family. To each cycle one can associate a vector c in 
IRl^l by 

{0 if the edge e does not belong to the cycle, 
1 if the edge e belongs to the cycle and is positively oriented, 
— 1 if the edge e belongs to the cycle and is not positively oriented. 

A family of elementary cycles is independent if the corresponding vectors are linearly independent. 
The cyclomatic number of the graph is the maximal number of independent elementary cycles 
of the graph; it is independent of the orientation. 
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clear once we use the fact that the cyclomatic number of a planar graph (V, E) 
also corresponds to the number of bounded connected components of delimited 



by the edges of the graph (which are called finite faces in ||Be|| ; see, for example, 
Theorem 2 therein). Then ( p. 311 ) amounts to saying that to each finite cluster of n 
corresponds one and only one such finite component of n, which is straightforward 
to prove. 

Below, we will use an extension of these relations in the case of infinite graphs 
where n is a configuration satisfying the +-boundary conditions. To make 
sense of the above formulae, we will apply them to the restriction of this graph to 
the graph {V,B), where 1/ := {t G : rfi(t,A) < 1} and B = B{V). This will 
give us the relation we require, up to some constant independent of n. 

Remark: We have only considered simply connected sets A. Let us make a com- 
ment in the case of non-simply connected A. To be specific, suppose the set A is a 
square with some hole in the middle, and that we have (+, +)-boundary conditions 
in the spin model. Then the associated random-cluster measure can be defined 
similarly as was done above, but the corresponding +-boundary condition are such 
that the two disjoint components of the the set {b : b ^ B^{A)} must be connected 
by an extra open bond. This makes the graph non planar and the relation ( p. 31 



does not hold anymore. Thus, as was the case for the duality of the spin system, 
the condition that A is simply connected is essential. More general settings can be 



studied using techniques of algebraic topology as in LMeR . 



3.4 Duality in the random— cluster representation 

Let A be a bounded simply connected subset of Z^. 

Let n = (no-,Ilr) be a configuration of a- and r-bonds. The dual configuration is 
defined as 

= 1 - n^,b ■ (3.32) 

We emphasize the fact that we have exchanged a and r bonds, this is done for later 
convenience. 

Note that if n satisfies the (+, +)- (resp. free) boundary condition on A, then n 
satisfies the free (resp (+, +)-) boundary condition on A*, and that S+(A) = B{A*). 
Indeed, let's consider what happens to a single site of A during the process of going 
to the random-cluster representation and then to its dual. 



O 



Figure 1: Lattice transformations (only one site shown) 

We start with the Ashkin-Teller model in A with (+, +)-boundary condition. Let 
us consider some site t G A. We then define the corresponding GRC model with 
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(+, +)-boundary condition: here the bonds whose value is not fixed are all bonds 
whose boundary contains at least one site of A (among them, there are in particular 
the four bonds with an endpoint at t). Then we take its dual model with free 
boundary condition: now, the bonds whose value is not fixed are the dual bonds 
crossing the previous ones. But the four dual bonds around t are the boundary of a 
plaquette p* having t as its center (see fig|5). Hence, to each t G A, this associates a 
plaquette p*{t) having t as its center: this is just the definition we have given of A*. 
Therefore these four bonds belong to B{A*). Doing this for all sites of A, we obtain 
all oiB{A*). 

Using the preceding geometrical results, we can write 

+A 



E /(i^)- n «o(&) n ^Ab) n n x 

.^^ bG6+(A): bGB+(A): 6gB+(A): beB+(A): 

n6=(0,0) n6=(l,0) ni,=(0,l) nt={l,l) 

c E /fe) • ( n «o(&) n n n «-(^)) >< 



free b.c. in A* 



fee/3(A*): 
«6={1,1) 



beB(A*): 
nj=(l,0) 



beB(A*): 
nj,=(0,l) 



beZ?(A*): 
"6={0,0) 



^qNr{h\h)+\h,\^Na{h\K)+\fk 

C 



E fin) 

fi: 

free b.c. in A* 



«g={l.l) 



beB(A*): 
ng={l,0) 



6eB{A*): 
ng=(0,l) 



feGB(A*): 
ng=(0,0) 



n «-(^) n «-(^) n n «o(fe) 

(3.33) 



where C denote some constant independent on n, \n\ is the number of open bonds 
in h and we have introduced 



ao(fe) 


= C'a„r{b) , 


aa{b) 


= qrC'a„{b) , 


ar{b) 


= q„C'ar{b) , 


dar{b) 


= qaqrC'ao{b) , 






Qt 


= qa- 



C is a normalization constant such that ao + cio- + + ^io- 



1: 



(3.34) 



(3.35) 



Remark: 1) Note that the positivity of theses dual probabilities is a consequence of 
positivity of the initial probabilities. 

2) There are no other way to distribute the factors q^ and qr in ( |3.33|) . 



As a consequence we have 
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Proposition 3.2 Let A be a finite, simply connected subset of ^ . Let clq, a^, a-, 
and Gar be defined by \3.3^ . Then, for all f e J-'^, 



u^{f\ao,a„,ar,a„r,qa,qT) =i^i*{f\ao,a„,ar,aaT,qa,qT), (3.36) 
where f{n) := f{n{n)). 

A natural question to ask is: do the transformations between the spin and random- 
cluster models and the duality transformations commute in the Ashkin-Teller case? 
The answer is yes, as is shown below. 



3.5 Commutativity of the dualities and RC transformation 

We want to compare the dual of the generalized random-cluster model associated to 
the Ashkin-Teller model with (+, +)-boundary conditions with the random-cluster 
model associated to the dual of this Ashkin-Teller model. 



For the dual of the random-cluster model, we obtain (see (p.34| ) and ( p.l4| )) 

1 -jarij^+jr) + jcrjr 



ao 



^ + jariJa + jr) + jajr 
^ + jarUa + jr) + jajr 



^ + jarUa + jr) + jajr ' 

So-T = 1 — cio — So- ~ Q^r ; (3.37) 
where j^ = e~'^^'' , jr = e"^"^^, j„r = e''^'^"'" ■ 

For the random-cluster of the dual model, we find (see (|2.11|) and ( p.l4| )) 

, _ l + st 

~ l + stV 
, _ (t-/)(l-s) 
l + stl ' 

{s-m-t) 



l + stl 

= I -a* -a: -a*. (3.38) 



where I, s, t have been defined in (^.6]). Using the relations: 



1 -ja ^ _ ^- jr jo 



Jar 



Jar 



^+ja l+jr 1 + jo 

it is easy to see that the quantities defined in (|3.37|) and ( p.38|) are in fact the same. 



We have already checked that both resulting models are defined on the same lattice 
(see section p.4[), so we can conclude that the following diagram is commutative: 
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AT 
RC 



AT* 

TIC 
RC* 



Here TIC denotes the random-cluster transformation and * the duahties. 



We now turn to the properties of the generahzed random-cluster model. 



4 Properties of the GRC measure 
4.1 FKG inequalities 

In this section we are going to show that the measures of the generalized random- 
cluster model which have been introduced are FKG. No hypothesis on A, except its 
boundedness, is required. 

We partition the bonds into two classes: 

B> = {6 : a^rib)ao{b) > a^{b)ar{b)} , 

B< = {6 : a,r{b)ao{b) < a,(6)a,(6)} . (4.1) 
We introduce the following partial order on{0,l}x{0,l}: 

(0, 0) ^ (0, 1) ^ (1, 1), (0, 0) ^ (1, 0) ^ (1, 1) , (4.2) 
for bonds in B>, and 

(0, 1) ^ (1, 1) ^ (1, 0), (0, 1) ^ (0, 0) ^ (1, 0) , (4.3) 

for bonds in B<. 

For the generahzed random-cluster associated to the Ashkin-Teller model, it is easy 
to see that all bonds will be in B> if Jar > and in B< if Jar < 0. 

Definition 4.1 Letm andn he two configurations, m is said to dominate n, myn, 
if mi, h nb, Mb. 

Definition 4.2 A function f is said to he increasing if m>z n=^ fiuk) > fin)- 1^ 
is said to he decreasing if —f is increasing. 

Example: N„{n\ts) and N^{n\K) + X]feGB< ^T,b are decreasing functions for the order 
defined above, while N^{n\A.) + X];,no-,6 and N-,-{n\A.) + Z^bes^ J^T.fc are increasing 
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functions. 

Let us just look at the case Nr{n\A) + J2beB> ^T,b- It is sufficient to consider two 
configurations nyfi, differing only by one r-bond at b. 

There are two cases: either 6 G i3>, or 6 G i3<. In the first case the r-bond is missing 
in n and hence we have 

iV.(n|A)+ = Nr{n\A)+ ^ n.,6+(iV,(n|A)-iV,(n|A))-l < N,{n\A)+ ^ , 

b€B> b€B> b&B> 

since |A^r(^|A) - Nr{n\A)\ < 1. 

If 6 G -B<, then the bond is missing in n and 

Nr{n\^)+ J2 <b = Nrin\A)+ J2 rir,6+(A^r(^|A)-iV,(n|A)) < iV,(n|A)+ ^ n,,6 . 

Indeed, the r-bond links two sites already in the same r-cluster and therefore 
the number of such clusters doesn't change, or it links two different clusters and 
A^^(n|A) - A^^(n|A) = -1. 



Definition 4.3 A measure is said to be FKG if fJ^ifg) > fi{f)fi{g), for all in- 
creasing functions f and g. 



Lemma 4.1 The generalized percolation measure Xb is FKG for the partial order 
introduced above. 



Proof. It is sufficient to check that (see |[FKG|| ) 



AB(nVn')AB(n An') > \b{ii)\b{j]!.) 

where a V fe denotes the least upper bound of a and 6, while a A 6 denotes their 
greatest lower bound. As is a product-measure, it is sufficient to verify this for 
each bond, which is straightforward. The only non-trivial inequalities are: 

\((1,1))A,((0,0))>A,((1,0))\((0,1)), 
for bonds in i3>, but this is satisfied by definition of this class of bonds; and 

A,((l,0))A,((0,l))>Afe((l,l))Afe((0,0)), 

for bonds in i5<, which is also true by definition. □ 



Proposition 4.1 Suppose ^o- > 1 ond gr > 1- Then the random-cluster measure is 
FKG for the partial order introduced above. 



Proof. It is sufficient to check (see |FKG[) that 



Af^(nVn'|A)+Ara{nAn'|A) iVr(nVn'|A)+Afr(nAn'|A) „Na{n\K)+Na(n;\A) Nr(n\h)+Nr{n!\^) 
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which can be proved exactly as in the case of the usual random-cluster model 
||ACCN|| . The only thing to observe is that 

A^^(nVn'|A) = N„{n9n^\A) , 
N^{nA]i\A) = N„{nAn\A), 
Nr{nAn\A) = Nr{n9n^\A) , 
Nr{n\/n\A) = Nr{n/\n!\A) , 

(4.4) 

where V and A denote the order induced by setting the order ( [4.2| ) at all bonds. 
Hence, 

A^<t(21V?i'|A) + A/'^(nAn'|A) = A/'^(nVn'|A) + A^^(nAn'|A) , 
Nr{nyn^\A) + Nr{nAn!\A) = A^^(nVn'|A) + A^^(nAn'|A) . 

n 

As direct elementary applications of these inequalities to the Ashkin-Teller model, 
we have 

|^l{(T^c^j) = z/;(i A j) > A and j A A^) 

> uli^^A^)uXiJ^A^)=^,lia,)^^lia,), (4.5) 

which is nothing more than one of Griffith's inequalities. Here o means any boundary 
conditions for the Ashkin-Teller model and • the corresponding boundary conditions 
for the random-cluster model. 

More interesting is the following inequality, which holds in the case J„r < (for 
which we cannot use Griffith's inequalities), 

(4.6) 

where we have used the fact that k^l is decreasing for the order ^. 
More generally, we have, for negative J„r, 

f^A{(rA(yB) > /iA(0"A)/iA(c^B) , 
f^AicTATB) < fiAicTA)f^A{TB) , 

as can be easily verified (this is true for free, as well as for (+, +)-boundary condi- 
tions) . 

4.2 Comparison inequalities 

There is a class of inequalities in the usual random-cluster model which is very 
interesting: they allow one to compare the probability of an event for different 
values of q and of the probability of occupation. It is possible to generalize these 
inequalities here, as shown now. 
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We consider two random-cluster measures, iy\ and z/a, with parameters flg, CLfjy CL-J-, 
(^cTT, 9(T, Qt and cio, (la, clt, o,aT, (ja, Qr , respectivcly. What is the relation between 

the probabilities of monotonous events computed with these two measures? We will 
consider only two cases, but others can be proved in the same way. 

We introduce the following notations 

Pcy = pT = —, (4.7) 
ga Qr 

Oq da dj- OifjT / . „\ 

"0 = — , Oct = — , OLt = — , "(TT = • (4.8) 

QiQ Qju Qjf 

We can now formulate our first inequality. 

Lemma 4.2 Suppose q^, qr > 1, and 

qa <qa, qr <qT, (4.9) 

tto-T > max{aa,ar) > mm{aa,ar) > a;o,V6 G i3> , (4-10) 

Pt-Oo- > max(a;o-r, Pr^o) > min(Q;o-^, PrO^o) > c^t , E , (4.11) 

then, for every increasing function A, 



the monotonicity of this function under the above hypotheses. 



Proof. Let xiR) = We are going to write x in such a way as to make explicit 



> 



X 



where C > is a constant independent on the configuration. 

Now, using ( [4.9| ) and the fact that Na{n\A) and A^T-(n|A) + X]f,GB< ^r,b are decreasing, 
it is easy to see that x "wiH be increasing if what is inside the brackets is increasing; 
and this will be true if it is true for each bond. This can be easily checked. We just 
consider two examples, since the other cases can be treated in the same way. 
Let us first verify that the expression in the first brackets is not decreasing when ni, 
increases from (1,0) to (1,1). In the first case this expression equals a^/aa, while 
in the second it equals aar/O'crT, which is not smaller by hypothesis. 
Let's now show that the expression in the second brackets does not decrease when 
Ufy increases from (1, 1) to (1, 0). But this amounts to aar/o-ar < qrO'a/qTO'a which is 
true by hypothesis. 

Doing the same computation for the other cases, we finally obtain 





z/a(A) = u^{A\x) > z?a(A), 
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by FKG and the fact that Pa and ua are normahzed. 

n 

We now give a second inequahty, 
Lemma 4.3 Suppose q^, Qt > 1; aiT'd 

qa>qa, qT>qT, (4-12) 

«o-T > max(pT-ao-, Pct^t) > niin(p^ao-, Po-ctr) > PaPrC^o , E , (4.13) 
> max(a^^,p^ao) > min(a(^^, p<^ao) > Paar,\/h e B< , (4.14) 

then, for every increasing function A, 

M^) > M^) ■ 

Proof. As above, using the fact that A''o-(?i|A) + \n^\ and Nr{n\A) + J2b€B> '^r,b are 
increasing. □ 

Remark: As said before, other such inequahties can be proved in exactly the same 
way, for example when g^r increases but q^ decreases. 

As a simple application of these inequalities, we prove inequalities relating the gen- 
eralized random-cluster model to the usual one. 

Lemma 4.4 Suppose q^, qr ^ ^, f & ^a, increasing, then 

PA{fa\Pl, qa) < MflQ^^ Qr) < pA{fa\p2, qa) , 

where 

qr{ao + a^) +ar + a 

The same kind of relations holds for f E Tr- 
Proof. 

^k{f\(la, Qt) < i^AiflQa, 1) = PkUAp = + a^T, qa) , 
where we used lemma |4.2| and lemma p.2|. In a similar way. 



MflQ^^ Qr) > MflQ^^ Qr = 1, ao = = Or = ^, a^r = 

= PAifaliqraa + a^r) / {qriO-O + Cia) + Or + aar),?^) , 

where N = qr{ao + a^j) + a-r + a^rr is the normalization of the new probabilities, and 



we used lemma and lemma 3.2. □ 



From this lemma we can for example obtain estimates on the location of the critical 
lines of the Ashkin-Teller model in the sector J^j. > Jcr = =: ^ > 0, in which the 



Random-cluster representation of Ashkin-Teller model 



24 




Figure 2: Schematic representation of the phase diagram in the plane tanh Jo-r versus 
tanh = tanh J^. We have indicated the critical lines (solid lines), and the self- 
dual line (which coincides with the solid line up to the splitting and then follows the 
dashed line). The shaded region corresponds to the set of parameters at which the 
random-cluster representation is not available. The estimates of the location of the 
critical line are also shown (dotted lines). 



critical line splits into two parts, the first one corresponding to the ordering of the 
ai and r^, and the other one to the ordering of the 9i = aiTi (see fig. ||). Let's study 
this second line. 

We first make the change of variables (cr, r) i-^ (^jt)- The previous lemma then 
gives, specializing to the increasing event {i^^A^}, 

PA(i^A1(2ae + aer)/(l + ao + ae),2) < z/a(z ^ A'=|2, 2) < pA(i ^ A>e + a^r, 2) , 
which can be rewritten in terms of the original Ashkin-Teller and Ising measures: 

r/A(^i| Jl) < ^iKie^i) < VAiOi\J2) , (4.15) 

where r]\{-\J) denotes the Ising measure with coupling constant J, and the coupling 
constants Ji and J2 are given by 

Jl = Jut + ^ lncosh(2J) , 

J2 = J + Jar ■ 

The upper bound in ( [4.15|) is easily obtained using GKS inequalities. This seems 
not to be the case for the lower bound. 

Since we know the critical temperature of these two Ising models, we obtain the 
following results (remember Oi = aiTi): 

Jar < Jfsing ' J ^ (^t) = , (4.16) 

Jar > Jfsing " 1 cosh(2 J) fx^ar) > , (4.17) 
where Jfsing = |argsinhl is Ising critical temperature. 
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Remark: Such a behaviour cannot occur in the sector Ja = Jt ^ Jar > 0. Indeed, 
in this case, 

/iA((Ti) = z/A(i A A^) > z/A(i A A'^ and i A A"^) = ^^{9i) 
which imphes fj,{9i) = -v^ = 0(= f^{Ti)). 



5 Possible extensions of the model 

The model above can be extended in several directions. 

The existence of the random-cluster representation and its properties (except dual- 
ity) do not depend on the particular structure of Z^. In fact all this can be shown 
to remain valid for an arbitrary finite subgraph A of some simple graph Q, applying 
exactly the same techniques as those used in this paper. 

The second direction in which the model may be extended is the following. Suppose 
we have possibly different Potts models on A, interacting through the following 
Hamiltonian 

N k 

^ = -E{E E 4"'-"^n(W*-i)}, (5-1) 

<ij> fc=l ri<...<rfc t=l 

which is an obvious generalization of ( ^.27] ) (af G {1, g^}). We can then define a 
generalized percolation measure {n := (n^, ...,nfc)) 

Mn):= n n (5-2) 

AC{1,...,N} bGB: 

where rij e {0, and 

N 

a^ = exp(E E - E «^ (5-3) 

k=l rl<...<r^.■. BCA 
r^^A,Mi B^A 

for all A C {!,..., AT}. 

These coefficients will then be positive under suitable conditions on the coupling 

constants, so that they can be interpreted as probabilities. 

It is then possible to define a generalized random-cluster model: 



if n satisfies the *-b.c. on A, 



^^in|gi,...,g7vj := <i E,AW)(n)nZr?^ 

otherwise. 

(5.4) 

where -k denotes boundary condition, and Nk{n\K) is the number of clusters of type 
k in the configuration n, i.e. A^fc(?i|A) := A^(n^|A). 

It will then be possible, introducing enough classes of bonds, to prove again FKG 
inequalities and then comparison inequalities. 

Again a proposition analogous to proposition |3.1| holds for these new models. 



Random-cluster representation of Ashkin-Teller model 



26 



6 Conclusion 

In this paper we have defined a generahzed random-cluster model and shown how it 
is related to the usual random-cluster model and to the Ashkin-Teller model. This 
new model still possesses the main properties of the usual random-cluster model, 
namely FKG inequalities, comparison inequalities and a duality transformation com- 
muting with the duality transformation of the Ashkin-Teller model. 
Only direct applications of the obtained inequalities have been given (correlation 
inequahties, inequalities relating the generahzed random-cluster model to the usual 
one, and estimates for the critical lines of the Ashkin-Teller model), however many 
known results about the random-cluster model can be extended in a straightforward 
way. One of our motivations was to develop tools which have been shown to be very 



useful in the study of large deviations in the Ising model (see e.g. [|, Pi]) 



7 Appendix 



Proof of Proposition |2.1J . The equality of the two partition functions follows 
from relations ( p.ll| ) and comparison of ( |2.4| ) and ( |2.1CI| ). Note that the summation 
is over all families of (compatible) closed contours without further constraints. This 
is the case because A is simply connected (see section p.l| ). 

1. We first show that ( 2.11|) is well defined, that is, that the functions S, T and L 
are strictly positive for any given triple {J*, J*, J*^) G V. 

This is obvious if J*^ > 0, so that we only consider the case J*^ < 0. In this case, 
we have 



s > st"^ > -It S > , 
t>ts^ > -sl^T >0, 
I > -st ^ L > 0. 



2. For every triple {J*, J*, J*^) G V, we can solve (|2.11| ) and get a unique triple 

( , Jfy J rrr ) • 

3. We now show that the map just defined on V 



i^Jfji '^Ti Jut) ' ^ (^CT) Jti Ja 



(7.1) 



takes its values in T>. 
3«ci J (J ^ J J- 



3.b X > 



3.C X > J„ 



J„>J^^ <l^T<S^s{l-l)>t{l-l) 



Jr>0^ SL<T ^ {1- f)ts^ < (1 - f)t . 



Jr > Jar ^ L <T ^ t{l - S) < /(I - s) . 
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3.d tanhJo-T > — tanhJo-tanhJ,- 

We use the following elementary result 

I - a 1-/31-7 

tanha > -tanli6 tanhc > ^ + 7) < 1 + ^7 , (7.2) 

1 + a l + pl + 7 

which holds for all triple of real numbers a, b and c, and a = e~^", f3 = e~^^, 7 = e~^'^. 
This gives 

tanhJ^^ > -tanhJ^tanhJ^ ^ e"^'^''" (e"^-^"^ + e"^-^") < 1 + e''^'^'^^^'^^^ 

^ S + T <1 + L 

^ /(l-s)(l-t) > -(1 -s)(l -t). 



4. We now prove that ( [7. 1|) is one-to-one. It is sufficient to show that for any triple 
( Jcr, Jt, Jar) € we can define a triple (s, t,l) G (] — 1, 1[)^ (see (p. 61)) and that the 
corresponding triple (J*, J*, J*^) G V. 
We claim that (s, t, /) is given by 

s = {I + - T^ - - [{I + - - LY - 4(5 - TL)2]5)/(2(^ - TL)) , 
t = (1 + - - - [(1 + - - LY - 4(T - 5L)2]5)/(2(r - SL)) , 
/ = (1 + - 52 - - [(1 + - 5' - T')^ - 4(L - 5T)2]^)/(2(L - ST)) . 

(7.3) 

4. a Let us verify that the quantities inside the square brackets are positive. 
4.a.l + T2)2 > 4(L - ST)^ 

We have 

tanhJ^^ > -tanhJ^tanhJ^ ^ ^~2Jar (^^-2J„ _^ e"^-'") < 1 + 6"^^'^'^+''") . 

^ S + T < 1 + L 
^ 5^ + - 1 - < 2(L - 5T) 



where we have used ( [7.2|) and the fact that S", T, L are positive. Now if < 
then 

L-ST = e-2(J-+Jr)^i _ e-4-/<.r) < g , 

hence, 

(1 + - ^2 _ ^2^2 > _ ^^^2 _ 

On the other hand, if Jo-r > 0, we have 

which is equivalent to 

1- L> S -T ^1 + - -T^ >2{L - ST), 

where we have used the fact that S* > T if > Jr- This last expression finally 
gives 

(1 + - - T^)^ > 4(L - ST)^ . 

4.a.2 (1 + T2 - 52 - L2)2 > 4(T - 5L)2 
We have 

1 + - 52 - = 1 + - 52 - + 2(T2 - L^) > , 
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because - L'^ = q-^J^ (^^-^J^r _ g-4J.^ > if > J^r- 
On the other hand, 

J^>O^T-SL>0. 

Thus, 

+ > 4(r - SLf ^{T + lf>{S-Lf^T+l>S-L, 

which holds if > J^. Then use I + L - S + T > S + T - S + T > 
4.a.3 (1 + 52 - T2 - L^f > A{S ~ TLf 

Again l + S^-T^ -L"^ ^1 + L^-S^-T^ + 2{S^ - L^) > 0, and 5 - TL > if 
> 0. So 

+ L^f > A{S - TLf ^ S + 1>T - L, 

using the fact that J^-r < Jt- The claim follows from S+l+L—T > S+T+S—T > 0. 
4.b We now prove that s, t, / e] — 1, 1[. 
4.b.l s > 

As 5" > LT (see 4.a.3.), it is enough to show that 

1 + 52 _ y2 _ ^2 _ p ^ ^2 _ ^2 _ ^2^)2 _ _ ^^^2|i > Q ^ 

but this is obvious. 
4.b.2 s < 1 

This is equivalent to show that 

4{S - LT)[2{S - LT) + L^)] < , 

which is a consequence of the above results (see 4.a.3.). 
4.b.3 l>t>0 

This is proved in the same way as for s. 
4.b.4 J^^ > ^ 1 > / > 

L — TS is positive and we obtain the same kind of relations as for s. 
4.b.5 J^^ < ^ > / > -1 

This time we have L — TS < 0, which gives the results in the same way as before. 
4.C It remains to show that {J*, J*, J*^) e V. 

We have already seen that J* > (see 4.b.l.), J* > (sec 4.b.3.), so we just have 
to prove that J* > J*, J* > J*^. and tanhJ*^ > -tanhJ*tanhJ*. 

4.C.1 j; > j; 



4.C.2 J* > J*^ 

In the same way. 



J*>J* ^ s>t 



« (i^Ml_ti)>o 

1 + stl 

^ S>T 
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4.C.3 tanhJ*^ > -tanhJ*tanhJ;! 

tanhJ*^ > -tanhJ*tanhj; ^ I > -st ^ L> 0. 

4.d The fact that (|7.3|) are solutions of (|2.7| ) can be checked by exphcit substitution. 

□ 
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